Suppression of Tc in superconducting amorphous wires 
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The suppression of the mean field temperature of the superconducting transition, Tc, in homo- 
geneous amorphous wires is studied. We develop a theory that gives Tc in situations when the 
dynamically enhanced Coulomb repulsion competes with the contact attraction. The theory accu- 
rately describes recent experiments on Tc-suppression in superconducting wires, after a procedure 
that minimizes the role of nonuniversal mechanisms infiuencing Tc is applied. 
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Disorder suppresses the superconductivity transition 
in morphologically homogeneous superconductors 0-0] 
because the diffusive character of the electron motion in 
dirty systems makes the Coulomb interaction more effec- 
tive 0. As a result, the attraction between the electrons 
in Cooper pairs becomes weaker, and the transition tem- 
perature, Tc, is lowered. In two dimensions {2D) the 
influence of disorder on T^ can be studied systematically 
by varying the film thickness d [|6|-g|. In uniform films 
Tc, being well defined, is suppressed as the sheet resis- 
tance, Rn, increases with decreasing d. (For a review 
see Ref. [^.) When the geometry of the sample is such 
that its dimension is lowered towards the one-dimensional 
(1_D) limit, the suppression of superconductivity should 
become more pronounced pO|] . 

Recently, efforts have been made |ll|-p^ to extend the 
experiments in films to narrow wires by fabricating a se- 
ries of amorphous Pb wires of different thicknesses and 
widths. It has been found that the Tc-suppression be- 
comes stronger as the wires' width reduces below lOOOA. 
The experiment of Refs. ||l2,n3[ is in the crossover region 
from 2D to ID. Actually, the wires are in the ID limit 
as far as superconducting fiuctuations are concerned p4[ , 
but they are in the crossover region from 2D to ID with 
respect to the diffusive motion of the electrons. 

From the theoretical point of view, the problem of Tc- 
suppression in ID wires is rather intriguing. As is well 
known, the superconductivity transition is determined by 
a series of logarithmically divergent terms describing the 
electron scattering in the two-particle Cooper channel. In 
2D systems the corrections due to the electron-electron 
(e-e) interactions combined with disorder are logarith- 
mically divergent as well g]. As the whole problem is 
controlled by logarithmic singularities, it can be stud- 
ied by renormalization group (RG) methods [Q. In ID, 
due to the reduced dimensionality, the effect of e-e in- 
teractions is more singular. It produces corrections that 
diverge as the square root of the frequency. The presence 
of two types of singularities demands a special analysis in 
the calculation of Tc- In this paper we develop a theory 
that describes adequately the effect of the dynamically 
enhanced e-e interaction on Tc in the crossover region 



from 2D to ID and perform a detailed comparison with 
the experiment. 

The mean field temperature, Tc, is defined as the tem- 
perature at which the electron scattering amplitude in 
the Cooper channel, Fc, becomes infinite. Fluctuations of 
the superconductivity order parameter lead to a broad- 
ening of the phase transition. However, its mean field 
temperature can be found experimentally by fitting the 
upper part of the resistive transition to the Aslamazov- 
Larkin theory [[l4| . The diagrammatic representation of 
the amplitude Fc is shown in Fig [l[ In addition to the 
contact BCS-interaction amplitude 7, the terms arising 
as a result of the interplay of the Coulomb interaction and 
disorder are also included in the Cooper ladder-diagram 
series. [The impurity scattering does not influence the 
e-e interaction mediated by phonons because in the long 
wavelength limit the lattice defects oscillate together with 
the ions |l6[.] The resulting equation for Fc is: 

rc(e„,e/) = -|7l+iA(e„-He/) 



-2t:TJ2 [-h\ + tA{en + em)] 



m=0 



X — Fc(e^ 



,£/), 



(1) 



where em = 2TrT{m + 1/2) is the Matsubara frequency, 
and the summation over m is limited by M ~ (2ttTt) 
In this equation 7, the bare value of the amplitude Fc, 
is rcscaled in such a way that the Debye frequency as a 
cut off energy is substituted by t~^, the inverse of the 
scattering time. Then, 7 = 1/ In (Tcor/1.14) , where Tco 
is the temperature of the superconducting transition in 
the bulk limit. The parameter t = {e^ /2n'^h)Rn charac- 
terizes the level of disorder in a sample, where Rn is the 
sheet resistance. The amplitude A describing the com- 
bined action of the e-e interaction and disorder is given 
by 

AK) = u^ > ^ ^ , , ^ , , (2) 
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where a and L are the width and the length of the wire, 
respectively. The parameter u describes the amplitude of 



the e-e interaction when the momentum q transferred by 
this interaction is not too small compared with the trans- 
ferred frequency w, namely when q ^ q^j = y^tu/D. (As 
was explained in Refs. MM, the most divergent contri- 
butions from the region q < q^ cancel each other out. In 
this region of small momenta, the e-e interaction depends 
only on the frequency, and therefore it can be gauged 
out.) Next, for amorphous Pb films the spin-orbit scat- 
tering time is expected to be only a few times longer 
than the elastic scattering time and therefore the part 
of the e-e interaction related to spin density fluctuations 
can be neglected. In that case, we may take u to be the 
value of the screened Coulomb interaction amplitude in 
the region of momenta q ^ quj, which gives u ~ 1/2. 
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FIG. qI The diagrammatic equation for the scattering ampli- 
tude Tc in the Cooper channel. The block 7 denotes the 
BCS-interaction amplitude. The block tA describes the inter- 
play of the Coulomb interaction with disorder that leads to 
the suppression of Tc- The wavy line is the screened Coulomb 
interaction, dashed lines describe impurity scattering. 

In 2D the summation in Eq. (0) yields A(w„) ~ 
Mln(l/a;jjT). Therefore, Eq. (|l|) combines the usual BCS 
logarithms together with the ones arising due to disor- 
der. Unlike the ladder diagrams in the BCS-theory, the 
integrations in the different blocks of the diagrams in 
Fig. |l| cannot be factorized, because A (e„ -f em) matches 
the frequency arguments of two neighboring blocks. In 
order to solve this parquet-like equation with a logarith- 
mic accuracy one uses the approximation ln(z -f- z') = 
ln(max{z, z'}), see e.g. Ref. [|7|. Then, it is possible to 
apply the "maximum section" method. This procedure 
leads to the RG equation for the amplitude Tc{e, e) P,|l5|: 
dTc/dlg = ut — r^, where l^ = ln(l/eT). The integration 
of the RG equation gives the suppression of Tc by the 
Coulomb interaction in 2D disordered systems: 



TcoJ \l\ 



1 , 1 + V^/|7| 



In 
2Vut 1 



ut/ I7I 



(3) 



This formula accurately describes the experimental re- 
sults in MoGe films H, with u = 1/2 and using only one 
fitting parameter, 7 H. 

In ID the result of the summation in Eq. (ph yields a 



square root singularity in the amplitude A(ti;„). When 
one deals with singularities stronger than logarith- 
mic ones, the approximations of the maximum section 
method cease to be valid, and a different method should 
be invented. In this Letter we treat the problem of find- 
ing Tc from Eq. (nl) as a sort of an eigenvalue problem, 
which leads to an implicit equation for Tc- To see this, 
we will consider rc(en,em) as the matrix elements of a 
matrix Tc, and will write the solution of Eq. (m for Tc in 
matrix notations: 



r. 



;^(/-|7|n) 'r^(-|7|i + tA) 



(4) 



\-l^ 



4A]e ^'^,J„m = Snrnin + 

>~nm = 1, and / is a unit matrix. 



Here n(T) = e-^/'^[\ - |7| 
1/2), A„„ = A(e„-fe„), 1, 

Eq. (0) is written in such a form that 11 is a symmetric 
matrix. Notice, that the dependence of II on the tem- 
perature T is not only through the dependence of A on 
the Matsubara frequencies, but also through the matrix 
rank M — {2'kTt)~ . The amplitude F^ diverges when 
the temperature is such that one of the eigenvalues of 
the matrix II(r) is equal to |7|^^, i.e., at T = Tc the 
equation 



|7r^/-ii(rc) |*) = o 



(5) 



holds. Thus, the equation determining Tc can be ob- 
tained from an eigenvalue problem. [One can also obtain 
an equation for Tc by considering a BSC-like gap equation 
with frequency dependent interaction vertex, — I7I -l-tA.] 



The matrix elements of 11 = 11" 



n^ 



ni. 



-t[{n 



,„ = [(n + l/2)(m + l/2)]-l/^ 

+ l/2)(m+l/2)]-i/2|^|-iA(e„ + e, 



(6) 



As the matrix elements 11°^ are factorized with respect 
to n and m, all the eigenvalues of the matrix 11°, ex- 
cept one, are degenerate and equal to zero. The eigen- 
vector corresponding to the nonzero eigenvalue is '^^ — 
c/i/n+ 1/2, and the equation |7|-i*° = E,„n°™*^ 
leads to the BCS relation for Tcq: 



M 



\l\-' = loiTco), loiT) ^ J2 



1 



m=0 



1/2 



In- 



1.14 
Y7' 



(7) 



Our strategy now will be to calculate the corrections to 
this eigenvalue perturbatively in 11^ (notice that 11^ oc t), 
and in this way to get an implicit equation for Tc- Since 
Ft is symmetric we can perform this program using a 
standard perturbation theory: 

\j\-^^l,{T) + hiT) + hiT) + ... (8) 

The first order term can be obtained straightforwardly 

;i = (*°|fiik°\ = --^S2(r), 



S2(r) 



^o|7l 
A (e„ + e„ 



M 

^ fn+l/2)(m+l/2)- 



(9) 



The prefactor l/lg appears in li because the normaliza- 
tion factor c of the eigenvector ^° is equal to I/^/Iq. 
Since all the eigenvalues of the operator 11'^ are degen- 
erate except the one under studying, it is also possible 
to find the higher order corrections using only the eigen- 
vector 1 4"°), without involving other eigenvectors. We 
demonstrate it here for the second order term, but a 
generalization to higher orders is straightforward. In the 
second order 
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where 



^0 n^n^ 



^°) = 



loh\ 



rS3(T), 



y (T) = V A(e» -F efc) A (efc + em) 



ifc=0 



{n + l/2){m+l/2){k + l/2)' 



(11) 



Inverting Eq. (jg) perturbatively in t and having in mind 
that |7|Zo(rco) = 1, we find 



In 



Tc 



-tME3(r,o) + r„ 



Tco 



dT 




(12) 



T=T, 



Since Eq. (12) gives an approximation for Iii{Tc/Tcq), 



while the measured quantity in experiments is Tc/T^q, the 
first two terms of the perturbative series are sufficient for 
the description of the Tc suppression, if the parameter t 
is not too close to a critical value where Tc vanishes. [The 
parameter t should be inside the radius of convergence of 
the series (12). Outside this radius the superconductiv- 
ity is completely suppressed.] In the 2D case Eq. ([l^ ) 
reproduces the first two terms of the expansion of the 
right hand side of Eq. (||) in powers of ui/7^: 



In 



Tco 



I 



;^i(2n-M)7 V7^ 



(13) 



We note that expansion (g_3|) does not contain a term 
oc i^/7^. There are several diagrams that gives contri- 
butions to that order, however, finally they cancel each 
other fia] . The main advantage of Eq. dl2) is that it is 
not restricted to a logarithmic accuracy, and can be ap- 
plied to the description of the crossover from 2D to ID 
systems. 

In the experiment of Xiong et al. ||l^ the mean field 
temperature of the superconducting transition, Tc, has 
been measured systematically for uniform Pb wires of 
various widths. The effective strength of the disorder 



characterized by i?n has been controlled by the wire 
thickness d. Before going to a detailed comparison of 
the theory with the experiment a few remarks are in or- 
der. The theory described above deals with the universal 
mechanism related to large scale distances that are of the 
order of the thermal length Lt oc ^J D jT). However, a 
number of other effects may also inffuence Tc when the 
thickness d is decreased. For example, the electron states 
quantization and the interaction of the electrons with the 
film's substrate can alter the parameters of the electron 
liquid. These nonuniversal effects of a short range origin 
are not addressed by the present theory. In some sys- 
tems, e.g., MoGe (see Ref. H,^), the discussed effect, 
originated from the interplay of the Coulomb interac- 
tion and disorder, is dominant, and the theoretical curve 
matches the experimental data at 2D. Unfortunately, as 
it is shown in Fig. g the theoretical curve for Ph films 
does not follow the experiment. This fact indicates that 
the effects of a short range physics are not negligible here. 
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FIG. y Comparison between the theory (solid line) and the 
experimental results ]1^ (dashed line) in Ph films. The fit- 
ting parameter 7 = —0.16 is determined from the initial slope 
of Tc{Ra)- The deviation between the theory and the exper- 
iment at large Rn shows that the interplay of the Coulomb 
interaction and disorder is not the only mechanism influencing 

To minimize the role of the nonuniversal effects, and 
make the comparison between the experiment and the- 
ory possible, we proceed in the following way. For each 
width the theoretical curve has been multiplied by the 
function x{Ra) = T^^ {Ra)ex/T^° {Ru)th- This func- 
tion is the ratio between the two curves presented in Fig. 
^. Here, the basic idea is that, because the widths of 
the wires are considerably larger than any microscopical 
scale, the influence of the short range effects on Tc in 
wires remains the same as in 2D films. In this way, we 
believe, the effect of the long range physics determining 
the crossover from 2D to ID systems can be captured by 
the present theory. To continue further we have to dis- 
cuss another complication. Unlike the case of 2D films, 
the limit of Rn -^ for a series of wires with a fixed 



width is somewhat ambiguous. For the discussed data, 
the extrapolation of Tc to the Umit Ru — > at a fixed 
width yields values that are not equal to the transition 
temperature in the bulk limit. (Moreover, the extrap- 
olated values behave in an irregular way as a function 
of the wire width.) Under this circumstance, we have 
normalized the theoretical curves in such a way that in 
the limit Ru — > the fitting curves for each width, a, 
start from the extrapolated Tco(a) — Tc{Ru -^ 0). After 
this normalization procedure and rescaling the theoreti- 
cal curves by x {Ru ) , the data for wires of different widths 
has been plotted together with the theoretical curves in 
Fig. ^ The fitting parameter 7 = —0.16, determined 
from the initial slope of the Tc{Ru) in 2D films, was the 
same for all wire widths. Notice that at Ru > 200057 
the suppression of Tc for the wire of the smallest width 
is about 1.5 times stronger than for the widest one. The 
agreement between theory, i.e. Eq. (Jl^), and experimen- 
tal data for all wires of different widths turned out to be 
very good. 



for the description of the crossover region between 2D 
and ID. By treating the problem as an eigenvalue prob- 
lem, we overcame the difficulties occurring because of 
the coexistence of different singularities in the equation 
determining Tc- In order to compare the available exper- 
imental results with the theory, we analyzed the data in 
a way that minimizes the role of nonuniversal effects of 
a short range origin. We believe that the theory could 
be tested further with superconducting wires fabricated 
from other materials, where the initial slope of Tc{Ru) is 
larger than in Pb films. 
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FIG. M Comparison between the theory (solid line) and the 
experimental data |l3| for wires of different width. The short 
range effects are excluded assuming that they are the same as 
in 2D. For each wire width, a, the theoretical curves and the 
experimental data points are normalized by the extrapolated 
value rco(a) = Tc{Ra -^ 0). For all width we use 7 = -0.16 
as in 2D. 

To summarize, we have developed a theory that de- 
scribes the suppression of the mean-field temperature of 
the superconducting transition in amorphous systems. 
The theory is based on the consideration of the suppres- 
sion of the contact attraction due to phonons, by the 
dynamically enhanced Coulomb repulsion. It is suitable 
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